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Abstract 



Cuajimalpa, 



Any surface is completely characterized by a metric and a symmetric tensor field satisfying a certain set of 
constraints, the Gauss-Codazzi-Mainardi (GCM) equations. These equations identify the latter tensor as the 
surface curvature. It should thus be possible, in principle, to phrase physical questions relating to a surface 
described by a Hamiltonian involving only surface degrees of freedom-be it a interface or a fluid membrane 
-in terms of these two tensors, without any explicit reference to the environment into which it is embedded: 
the surface itself is an emergent entity. Introducing Lagrange multipliers to impose the GCM equations as 
constraints on these variables, we demonstrate how the Euler-Lagrange equation describing the stationary states 
of any surface Hamiltonian is derived. The behavior of these multipliers is explored in detail for surfaces 
minimizing area. Using a soap film between two rings as an example, it is shown how singularities in these 
functions correlate directly with instabilities in equilibrium surfaces. This approach to interfaces and membranes 
has clear relevance to a number of problems in soft matter: in particular, this framework is ideally adapted to 
study the recently proposed programmed swelling of thin polymer sheets. 



^1 Introduction 

Surfaces occur as approximations of physical systems at almost all energy scales [lj. More often than one would 
^ej<pect the only relevant degrees of freedom are the ones associated with the geometric configuration of the surface 
T 7Ttfself and its behavior is described completely by a Hamiltonian or an action constructed using the geometrical 
invariants of this surface. This may be something as simple as the area-representing the energy of an interface or 
soap film [2j-or its relativistic analogue which represents the area of the worldsheet swept out in the course of 
Hrie evolution of a string, be it a fundamental extended object or-more conservatively-some effective description of 
one|l|4]. 

Typically the Hamiltonian defined on a surface, T : {w 1 ,?/ 2 } — > X(m 1 ,m 2 ), is constructed by forming suitable scalars 
using the induced metric g a b, the curvatures K a b and their covariant derivatives: 



H= / dAH[g ab , K, 
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for simplicity, we consider only surfaces embedded in three-dimensional Euclidean space, E 3 . The important point 
is that the functions X tend not to appear explicitly in H. Surface area with % = 1, depending on the metric 
through its determinant, dA = d 2 u^fg, is the simplest example. If the tensors g ab and K ab in Eq. ([!]) are to 
represent a surface, however, they will need to be consistent with the Gauss-Codazzi (GC) and Codazzi-Mainardi 
(CM) equations, 

C ± = 0, C a = 0, (2) 

where 

C± := H-K 2 + K ab K ab ; (3a) 
C a := V 6 (K ab - g ab K) , (3b) 

which occur as integrability conditions on the structure equations defining how the unit tangents and normals rotate 
as one moves along the surface. Here V a is the covariant derivative compatible with g ab ; 1Z is the corresponding Ricci 
scalar curvature and K represents the trace of K\, K = g ab K ab . Conversely, one of the cornerpieces of nineteenth 
century geometry is the assertion that any two tensor fields, g ab and K ab , satisfying Eqs. ([2j will represent some 
surface X, with induced metric g ab and extrinsic curvature K ab , unique up to Euclidean motions |5J. This will also 
be crucial. Indeed, even if H depended only on the metric, this metric knows there is an extrinsic curvature tagging 
along. 

In this paper, we will show that it is always possible to rephrase the variational properties of surfaces in terms 
of a theory of gravity involving a metric, coupled to an auxiliary field K ab , without any explicit reference to the 
embedding functions themselves: the surface itself is an emergent entity. In this framework Eqs. ^ are enforced 
by introducing Lagrange multipliers. This approach contrasts dramatically with the familiar approach in terms of 
harmonic maps, or its natural extension-when curvatures are involved-in terms of auxiliary variables jSlElE]: here, 
the surface does not materialize until these constraints are applied. 

A peculiarity of two-dimensional surfaces is that the multipliers assemble into a spatial vector field. If H depends 
only on the intrinsic geometry, this vector field can be identified in equilibrium as a generator of surface isometries; 
if H depends also on the curvature K ab , on the other hand, it is identified with a conformal transformation. The 
role of the multipliers themselves, however, is not to displace the surface. This identification is a two-dimensional 
accident: they represent the strength of the interaction coupling the tensor field K ab to the metric on the Riemannian 
manifold in the formation of the equilibrium surface. The surface Euler-Lagrange equations are derived by examining 
the flows generated by this vector field. Its behavior will be explored in detail for area minimizing surfaces. In the 
case of a catenoid bridging two rings, the relevant isometry will be identified explicitly, and the connection between 
the singularities in this vector field and the presence of instabilities emphasized. 

2 Look Ma, no X: Surface variational principles without surfaces 

Consider the following effective action or energy 

H c [g ab , K ab , A, A a ] = H[g ab , K ab ] + I[g ab , K ab , A, A a ] , (4) 

where 

I = \ f dAAC ± + X - J dAX a C a . (5) 
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The Lagrange multipliers fields A and A" enforce the GC and CM equations, Eqs. ([2]), as constraints on the variables 
g ab and K ab . In Eq. ^ one is now free to treat g ab and K ab as independent variables. The variation of H c is given 
by 

5H C = J dA (-^(T ab + T ab ) 5g ab + (H ab + U ab ) 5K ab + V a Q a ^ , (6) 

where the manifestly symmetric second rank tensors T ab and H ab , are associated with the variation of H with 
respect to g ab and K ab ; T ab and H ab are their counterparts for the constraining term /: 

T ab = -2^-, T ab = -2^- H ab = ^, u ab =^. (7) 
og a b og ab 6K ab 6K ab 

In Eq. ([6]) Q a represents all of the terms that have been collected in a divergence after integration by parts. 

The Euler-Lagrange equations for g ab and K ab describing the equilibrium states of the surface are given respectively 
by: 

T ab + T ab = Q. (g^ 

H ab + H ab = 0. (8b) 

supplemented with Eqs. (|2]). Eqs. (|8]) are analogues of the Einstein equations in general relativity. The technicalities 
of the variations with respect to g ab and K ab are themselves straightforward (see for example |8J). One identifies 

T ab = - l - (V a V 6 - g ab V 2 + TZ ab ) A + l - [V a (\ c K b c ) + V b {\ c K a c )] - l -V c [\ c K ab + g ab \ d K c d ] ; (9a) 
U ab = - (K ab - g ab K) A + - (V a X b + V b \ a ) - ^V c \ c g ab . (9b) 

The task is now to solve, if only implicitly, Eqs. ([8]) for the multipliers. This is facilitated by organizing the two 
tensors T ab and 1-L ab in a more geometrically transparent way. 

Introduce the Lie derivative along the vector field A a on the Riemannian manifold, which acts on the tensors g ab 
and K ab as follows 

Cxgab = VaAb + VbAa; (10a) 
C x K ab = (V c K ab + K ac V b + K bc V a )X c . (10b) 

and define analogues Ca for the scalar A 

C A g ab = 2K ab A; (11a) 
C A K ab = {-V a V b + K ac K c b )A. (lib) 

Let us cheat a moment by introducing the embedding functions X we claimed would play no role in the variational 
principle. Let e a = <9 a X be the tangent vectors adapted to the parametrization, and n the unit normal to the 
surface. One can then construct a space vector A = A a e a + An, with tangential components A a and normal 
component A. Now define 

^Adab = £a£W> + £\9ab ; (12a) 
C A K ab = C A K ab + C x K ab . (12b) 
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The induced metric g ab = e a ■ e b , and the extrinsic curvature tensor K ab = e a ■ V b n, then transform respectively by 



Eqs. (12) under the flow generated by the vector field A (see, for example, |9J). 



It should be stressed that neither of the definitions Eqs. ( |12| ) make any reference to the embedding functions, the 
identifications of g ab and K ab in terms of these functions, or the assembly of A and A a into a space vector; more 
importantly, despite the shorthand, it is not even appropriate to think of A as a spatial vector field. If this were 



the case the flow defined by Eqs. (12) would displace the surface geometry away from equilibrium. It is, in fact, a 
two-dimensional accident that a space vector can be constructed using the multipliers A a and A. In this context, 
the role played by A contrasts with the one played by the lapse and shifts in the Hamiltonian formulation of general 
relativity where the analogs of the GCM constraints for a spatial hypersurface embedded in a Riemannian manifold 
are the generators of normal and tangential deformations of this hypersurface [10]. Their role here is not to displace 
the surface: rather they are the generalized forces coupling the tensor fields g ab and K ab to form the induced metric 
and extrinsic curvature of the surface. 



The motivation for introducing Eqs. (12) is that it is now possible to cast the tensor T ab in the remarkably simple 
form 

T ab = \ (g ab K cd - g cd K ab ) C A g cd + l - (g ac g bd - g ah g cd ) C A K cd , (13) 
linear in C A g ab and C A K ab . In addition, H ab can be cast as 

n ab = \ {g ac g bd - g a b g cd ) C A g cd . (14) 

A useful identity: Using the intrinsic definition of 1Z, and its extrinsic counterpart implied by the GC equation, 
C± = 0, one identifies two equivalent expressions for C A 1Z: 

C A K = -2(n ab C A g ab + (K ab -Kg ab )C A K ab ), (15a) 

= ( W a \7 b -g ab V 2 -n ab )C A g ab . (15b) 



As a consequence, the projection of T , given by Eq. (13), on K ab can be cast completely in terms of C A g ab 



K ab T ab = i (g ab V 2 - V a V fe + K ac K b - g ab K cd K cd ) C A g ab . (16) 



The significance of this identity will soon be apparent. 



3 Two cases of interest: tension and bending 



Gravitational Impostors: 
in Eq. (§. 



Let us first consider a Hamiltonian depending only on the metric, so that "H = T-L[g a b] 



In this case H ab = so that Eq. (8b) implies that U 



lab 



0. 



The identity (14) in turn implies that the vector 

The 



0. 



A, treated as a space vector, can be identified as the generator of an isometry, in the sense that C A g cd 
identity ( 16 ) then implies that K ab T ab = 0. As an immediate consequence of Eq. (8a ), the Euler-Lagrange equation 



follows: a surprisingly short story once the role of A as generator of isometries is recognized. Notice 
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that T ab generally does not vanish; thus T ab ^ 0. Eq. (13) then implies that the isometry is non-trivial; for if 
^-\K a b 7^ 0, A cannot be a Euclidean motion. 

In particular, in the case % is some constant a, so that H is proportional to area, one identifies T ab = —ag ab , and 
the Euler-Lagrange equation reduces to K = 0. A familiar statement is recovered: the stationary states are minimal 
surfaces. 

Bending energy: A less simple example is provided by the Polyakov or Helfrich bending energy, quadratic in 
curvature, with H[g a b, Kab] = K 2 /2. It is the simplest non-topological conformal invariant of an embedded two- 
dimensional surface [llj; it also provides an extraordinarily robust mesoscopic description of fluid membranes [12J. 
Now T ab = K(2K ab - l/2g ab K) and H ab = Kg ab ; as a result of the latter, Eq. Q reads C A g ab = AKg ab . 
Thus A generates a conformal transformation, scaling locally with the mean curvature. Eq. (16) then implies that 
K ao T ab = (V 2 — K ao K ab )K so that the Euler-Lagrange shape equation is given by 

-V 2 + Tl- ^K 2 ^j K = 0, (17) 

a surprisingly pithy derivation that compares favorably with any of its better established counterparts, all the more 
so because this framework was not developed to compete on this level. 

The linearity of the Euler-Lagrange equations in the multipliers permits one to treat more complicated energies, the 
Helfrich Hamiltonian H = (K — K ) 2 /2 + a, with spontaneous curvature K and constrained area, for example. 



4 Multipliers for minimal surfaces 

It is curious that one never needed to identify the multiplier fields explicitly to isolate the surface Euler-Lagrange 
equations. If one were to stop here, however, would be a mistake: for in the role that they play in quantifying the 
forces necessary to constraint the tensor fields ([2]), the multipliers also signal when surface instabilities are present. 

In this section, the partial differential equations describing these fields will be determined. For simplicity, examine 
the area H = a with Euler-Lagrange equation K = 0. In general, the equation 

r a a = -^£A(vW • (is) 



follows by tracing over Eq.(13). Under the isometry A, (18) implies that the mean curvature changes by a constant: 
C\K = — 4cr. Combining this result with the contraction of Eq. (12b), £\K = (— V 2 + 1Z) A, one obtains 

(-V 2 + K) A = -4a . (19) 

The scalar A is determined independently of the vector field A a . The differential operator appearing here, Jz? = 
— V 2 + 1Z also makes an appearance in the second variation of area about any equilibrium geometry, which assumes 
the form 

S 2 A = [ cL4$^$, (20) 



where $ is the normal deformation of the surface. The existence of negative eigenvalues signals a mode of instability 
of the surface. 
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As discussed elsewhere [13J the appropriate boundary condition on A in Eq. (19) is A = 0. Its solution subject to 
this boundary condition is also unique. 



To complete the determination of A, note that the contraction of Eq. (12a) implies that V a A a = 0. The divergence 
of Eq. (12a) then reads 



(V 2 + l -K)\ a 



-2K nh V b A . 



(21) 



A sufficient boundary condition is A a = [15] . Given the function A, the solution of Eq. (21) is now unique. We 
will now show that the behavior of the multipliers correlate with the stability of the equilibrium surface. 

Example: Catenoid. We will examine the behavior of the multipliers on a catenoid bounded by two rings a 
fixed distance apart. Aligning the axis of symmetry along the Z axis, its radius and height, R(l) and Z(l), can be 
parameterizing in terms of arc-length I along its meridians (with / = on the neck of radius R , see Fig. [l{a)): 



R(l) = Vl + l 2 , Z(l) = arcsinhZ 



(22) 



where all lengths are measured in units of R . The principal curvatures along the parallels and meridians are 
Cm = — C± = 1/R 2 . or simplicity, consider a symmetric section of catenoid bounded by the parallel circles at 
I = ±L, with corresponding radius Rl and height ±Zl respectively. 



By symmetry A is axially symmetric. Eq. (19) then assumes the form 



(23) 



where the prime indicates a derivative with respect to arc-length, and A := A/(4ct_Rq). An exact solution of Eq. 
(23) exists. With the boundary conditions A(±L) = it is given by 

1 



A 



R 



\n(R + l) + l 



C n 



In (R + I) 



(24) 



with integration constant, 



Cn 



L ( LR L + \n(R L + L) 
4 \R L -L\n(R L + L) 



(25) 



Note that the global minimum, A , occurs at the neck where the curvature is highest. Furthermore A := A(0) = 
—Co diverges as L is increased to the value Lc = 1.50888 which occurs when Rl = L\yl(Rl + L) and the ratio of 
separation Kl = 2Zl to diameter of the rings Dl = IRl is h^jDi, = 0.66274. A is plotted as a function of I for 
several values of L in the interval [0,Lc] in Fig. [TJb). It is negative everywhere in this interval. 

The divergence of A at L = Lc correlates with the onset of an instability in the catenoid as a minimal surface 
(see Fig. [ljc)) . For let us expand A in terms of the eigenfunctions of the operator £6 ', A = ^ n C n Q n , where 
Jz?$ n = E n Q n , so that Eq. (19) reads ^ n S n C n $ n = — Aa. Let $ n be the normalized ground state with eigenvalue 
Eq. Then 



E C = -4cr / dA% 



(26) 



If $ is positive everywhere, the left hand side of Eq. (26) is manifestly negative. If L is small, the catenoid 
approximates a cylinder with positive E . This implies that Co is negative and thus so also is A, consistent with 
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the exact solution. As L — > L c , however, E — > 0. At this value of L Eq. (26) implies that C must diverge, 
so that A does also. Thus an unexpected bonus of this framework is a reformulation of the analysis of stability 
of minimal surfaces. L = Lc is the maximum value of the meridian length for which the catenoid is stable [14|. 
Beyond L = Lc, E becomes negative and it can be shown that A changes sign. 





(b) 




0.0 0.2 0.4 0.6 



(c) 



L/L c 



0.5 0.0 
I 

(d) 



Figure 1: (a) Catenoid between two identic rings, (b) A as a function of arc-length I for values of boundary arc-length 
L = 0.6,1,1.25. The gray scale increases with L. (c) Local extrema of the multiplier A/(4<ri?o) as a function of the 
boundary arc-length L. (d) A; as a function of arc-length I for the same values of L as in (a). 



Likewise, defining A Q 



A a /(4cri?Q), Eqs. (21) are given by 
1 



R 



2— A' 
R 2 



-A, 



0. 



(27a) 
(27b) 



Axial symmetry implies that the angular component vanishes:A^ = 0, consistent with the fact that the CM constraint 
equation along the polar direction vanishes identically: = 0. Thus, there is only a generalized force along the 
meridians. The corresponding component A; is plotted in Fig[T][d) for values of I in the interval [0, Lc]. It is an 
antisymmetric function of /, possesses two extrema, one maximum and one minimum, vanishing at the neck where 
I = 0. Like A, A^ diverges at the onset of instability at L = L C - 



5 Discussion 



In this paper it was shown how a surface can be treated as a Riemannian manifold endowed with a metric that 
couples to a symmetric tensor field. The GCM equations impose a constraint on these two fields. No direct reference 
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is made to the surface itself. 



This approach does not replace the more intuitive direct approaches, but nor was it intended to. Indeed, as we 
argued it has its shortcomings: it is not obvious, for example, how to address the Euclidean invariance of the surface 
Hamiltonian, never mind the conservation laws that it implies, when the surface and its background do not yet exist. 
The motivation, however, was never to compete with tried and tested methods: it was to establish a framework for 
studying surfaces that mimics gravity; the surface itself is an emergent entity. In the process, intriguing connections 
with a theory of metrics are revealed that are likely to be worth exploring. 

Relevant antecedents of this work can be found in Barbour, Foster and Murchadha's Relativity without relativity, 
Sorkin's treatment of field theory in Minkowski space, or Lomholt and Miao's discussion of the ambiguities associated 
with the GCM equations [15]. It also shares features with the framework, developed in [16J in the context of paper 
folding, for enforcing local geometrical constraints. In this context, the metric approach developed here-tweaked 
appropriately-is likely to be relevant to the fascinating current research on programmed swelling [17J. A toy model 
might involve minimizing bending energy, while constraining the metric g a b on the target geometry to coincide with 
an appropriate rescaling of the template metric g® b . 

This approach to interfaces and membranes has clear relevance to a number of problems in soft matter: fluctua- 
tions or membrane mediated interactions could be treated in a manner that sidesteps the difficulties of the height 
function representation. Numerically one could contemplate relaxing the GCM equations, but suppress violations in 
a controlled way by introducing large coupling constants. 

One curiosity and unexpected virtue of this framework is that the derivation of the surface Euler-Lagrange equation 
never requires the explicit determination of the Lagrange multipliers enforcing the constraints. These multipliers are, 
however, of considerable interest in their own right: it is they that quantify the strength of the coupling between 
the Riemannian metric and the symmetric tensor field shaping the manifold into a stationary state of the surface. A 
connection between conformal transformations and surface states has also emerged; its significance remains to be 
explored. More importantly, however, singularities in the multipliers correlate directly with instabilities in equilibrium 
surfaces. We have explored in some detail the behavior of these multipliers for surfaces minimizing area and, in 
particular, a soap film between two rings. 

Extending this framework to higher dimensional surfaces or non-trivial backgrounds is not entirely straightforward 
[13J. Unlike the two-dimensional case examined here, where the contractions of the GCM constraints completely 
encode their geometrical content, these constraints will need to be accommodated within the Hamiltonian in their full 
uncontracted glory. In particular, the fortuitous similarity with the ADM formulation of general relativity encountered 
here becomes an unreliable guidepost; the multiplier fields no longer assemble naturally into a vector field. What 
is more, the GCM equations will need to be supplemented with their Ricci counterpart if higher codimensions are 
contemplated [5J. 
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Appendix A Making connections and discovering shortfalls 



It is useful to compare this approach with a variational framework introduced by one of the authors several years ago 
which adopts a very different strategy |7J. In that approach H is again constructed using the metric and extrinsic 
curvature as independent variables. In contrast, however, these variables are connected to the embedding functions 
through the Gauss- Wei ngarten structure equations. One thus constructs the Hamiltonian 

H c = H[g ab , K ab ] + JdA [f • (e - d a X) + A^e, • n + A n (n 2 - 1)] 

+ J dA [A ab (K ab - e a • d b n) + X ab (g ab - e a • e 6 )] , (A.l) 

implementing the definitions of g ab and K ab in terms of the tangent vectors e a and the normal n, as well as the 
connection of the latter to X by introducing appropriate Lagrange multipliers. One is then free to treat each of 
these variables independently. In particular, the translational invariance of H implies the existence of a conserved 
stress tensor. In this framework, X only appears in the tangential constraint so that 

^H c = j dAVJ a ■ <5X (A.2) 

modulo a boundary term. Thus, in equilibrium, V a f a = 0, or the stress f a is conserved. f a is constructed using the 
remaining Euler-Lagrange equations. One finds [7], f a = f ab e b + f a n, where 

jab = T ab _ H ac R b ? ja = _y ^ab ? ^ 3) 

and T ab and H ab were defined in Eq. ([7|. It depends only on the geometry. This existence of this conservation law 
is not evident in the new framework. 

The normal projection of the conservation law reads 

V a f a -K ab f ab = 0; (A.4) 

its tangential counterparts V a f ab + K b a f a = 0, are the statement of reparametrization invariance. Notice that if 



"H depends only on g ab , Eqs. (A.4) reduces to the statement that —K ab T = 0. This also justifies the strategy 



that was adopted to identify the surface Euler-Lagrange equation. 
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